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ABSTRACT 

In the previous paper of this series Deshpande & Rankin (2001) reported results re- 
garding subpulse-drift phenomenon in pulsar B0943+10 at 430 MHz and 111 MHz. 
This study has led to the identification of a stable system of subbeams circulating 
around the magnetic axis of this star. Here, we present a single-pulse analysis of our 
observations of this pulsar at 35 MHz. The fluctuation properties seen at this low fre- 
quency, as well as our independent estimates of the number of subbeams required and 
their circulation time, agree remarkably well with the reported behavior at higher fre- 
quencies. We use the 'cartographic'-transform mapping technique developed in Paper-I 
to study the emission pattern in the polar region of this pulsar. The significance of 
our results in the context of radio emission mechanisms is also discussed. 

Key words: pulsars: general — pulsars: individual: B0943+10 — radiation mecha- 
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1 INTRODUCTION 

The foregoing paper by Deshpande & Rankin (2000; here- 
after Paper-I) reported detailed studies of the steady sub- 
pulse drift of pulsar B0943+10, using Arecibo polarimet- 
ric observations at 430 and 111 MHz. They introduced sev- 
eral new techniques of fluctuation-spectral analysis and ad- 
dressed the possible aliasing of the various features in an 
effort to determine their true frequencies. The paper further 
introduced a 'cartographic' transformation (and its inverse) 
which relates the observed pulse sequence to a rotating pat- 
tern of emission above the pulsar polar cap. This approach 
provided a direct means of determining some of the key pa- 
rameters characterizing this underlying source pattern and 
its stability — thus facilitating a fuller comparison with the 
predictions of pulsar emission theories (see also Deshpande 
& Rankin, 1999, hereafter DRa). 

While the Arecibo observations in Paper-I document 
a great deal about 0943+10's remarkable 'drift' behavior, 
they also demonstrate the need for further studies of differ- 
ent types and at various wavelengths. The 430-MHz obser- 
vations, for instance, sample only the outer radial periphery 
of the subbeams defining the overall hollow emission cone 
of the pulsar — and this is apparently a major reason for 
its unusually steep spectrum above 100 MHz. Observations 
at decameter wavelengths are particularly interesting, both 
because the pulsar is expected to be relatively brighter and 
also because they provide a means of exploring the char- 



acter of the larger low frequency conal beam size [often at- 
tributed to 'radius-to-frequency mapping' or 'RFM'; e.g. see 
Cordes, 1978, and the generally assumed dipolar field geom- 
etry]. At 35 MHz, then, we would expect our sightline to 
sample the region of emission more centrally, with the result 
that, the full radial extent of the subbeams would be 'seen' 
and the high frequency 'single' profile would be well-resolved 
into a 'double' form. The character of the low frequency 
emission pattern would also be of interest — both in its own 
right and as compared to those at higher frequencies — for 
what insights it might provide into the manner in which 
emission sites at different heights/frequencies are connected. 

Below, we report the results of our analysis of pulse- 
sequence observations of B0943+10 at 35 MHz. Our obser- 
vations and analysis procedures are outlined in the next 
section. Subsequent sections then describe the results of 
fluctuation-spectral analysis, the implied subbeam-system 
topology, and the derived polar emission maps [some pre- 
liminary results were reported earlier in Asgekar & Desh- 
pande (2000)]. 



2 OBSERVATIONS AND ANALYSIS 

Observations of this pulsar were carried out as a part of a 
new programme of pulsar observations (Asgekar & Desh- 
pande, 1999) initiated at the Gauribidanur Radio Observa- 
tory (Deshpande, Shevgaonkar & Shastry, 1989) with a new 
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data acquisition system (Deshpande, Ramkumar & Chan- 
drasekaran, 2000). The pulsar was observed for a typical 
duration of <; 1000 s in each session. Several such observa- 
tions were made during the spring of 1999. We discuss two 
data sets below, denoted as 'C and 'G', which were observed 
on February 7 and 13, respectively. 

From the signal voltages (over a 1-MHz bandwidth) 
sampled at the Nyquist rate with 2-bit, 4-level quantization, 
a time-series matrix of intensity with a desired spectral and 
temporal resolution (equivalent to data from a spectrome- 
ter) was obtained. After checking for, and excising, interfer- 
ence in the spectral as well as the time domains, the differ- 
ent spectral channel data were combined after appropriate 
dispersion-delay correction. The resulting time sequence was 
resampled at the desired intervals of the rotation phase of 
the pulsar to obtain a gated pulse sequence for a specified 
gate- width (typically ~ 100°). 

Our techniques of spectral analysis for the 35-MHz ob- 
servations are virtually identical to those used in Paper-I, so 
we will rely on its fuller descriptions of these methods and 
request that the reader refer to them there. 



Fluctuation Spectra and Aliasing 

Our 35- MHz observations lack the sensitivity to detect a ma- 
jority of the individual pulses in a sequence, let alone to iden- 
tify their drift pattern. Any such drift can best be investi- 
gated using longitude-resolved fluctuation (hereafter, 'LRF') 
spectra (Backer, 1973). Figure 1 then gives such a spec- 
trum for the G sequence. Note in particular the prominent 
peaks at 0.459 c/Pi (feature 1) and 0.027 c/Pi (feature 2). 
Feature 1 is associated with the observed subpulse drift at 
meter wavelengths; it is barely resolved here, implying a Q 
of >500. This clearly indicates that the underlying modula- 
tion process is also quite steady at decameter wavelengths. 
The modulation frequency is slightly different for the C ob- 
servation, but small variations (< 0.5%) were also noted by 
Deshpande & Rankin (Paper-I), so are not surprising. The 
important questions we wish to address are the effect of a 
possible aliasing of feature 1 and its relation to feature 2. 

To assess the issue of aliasing, we compute a harmonic- 
resolved fluctuation spectrum (hereafter, 'HRF' spectrum; 
see Paper-I: Figure 4) by Fourier transforming the continu- 
ous time series. Such a time series is reconstructed from the 
(gated) pulse sequence by setting the intensity in the off- 
pulse region to zero. A spectrum constructed in this manner 
for the G observation is displayed in Figure 2. Note that 
amplitude modulation would produce a pair of features in 
the HRF spectrum, symmetric about frequency 0.5 c/Pi and 
with nearly equal amplitudes; whereas, phase modulation 
would produce a strong asymmetry in the pair amplitudes, 
as exemplified in Figure 2. We can then see that the fluctu- 
ation feature at 0.541 c/Pi is associated with a phase modu- 
lation, and further appears as the first-order aliased feature 
1 at 0.459 c/Pi in the LRF spectrum of Figure 1. A similar 
conclusion regarding the nature of the modulation can be 
drawn by looking at the rate of change of the modulation 
phase with longitude (see Figure 5), which is 26° /° — too 
large for a mere amplitude modulation. 

We now fold the pulse sequence at the modulation pe- 
riod P3 in Figure 4. The drift-bands related to the two com- 
ponents are clearly visible and exhibit a phase difference 
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Figure 1. LRF spectrum of the 34.5-MHz G sequence, using 512 
pulses. The body of the figure gives the amplitude of the features 
in (arbitrary) grey-scale displayed on the right. The average pro- 
file of the pulsar is shown in the left-hand panel, and the integral 
spectrum is given at the bottom of the figure (see Papcr-I Figure 1 
for details). 
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Figure 2. HRF spectrum of the 34.5-MHz G sequence. The am- 
plitude of the frequency components at the fundamental rotation 
frequency (1/Pi) and its harmonics are plotted in the left-hand 
panel of the diagram. The body of the figure gives the amplitude 
of the other frequency components in the spectrum upto 40/Pi in 
(arbitrary) grey-scale, displayed on the right. The bottom panel 
shows the sum of the frequency components, collapsed onto a 
1/Pl interval. See Paper-I: Figure 4 for details. 
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Figure 3. HRF spectrum of the 34.5-MHz C sequence using 
128 pulses. The feature 2 (seen in G sequence) is not detectable 
here. 

between the modulation peaks. The 'bridge' region between 
the two average-profile components exhibits little discernible 
fluctuation power, so we could be missing one or more drift 
bands in this longitude interval. We therefore express the 
modulation phase difference between the two components 
as, A9 = 2n (0.8 + m) (where, m = 0,1,2..). If the modulation 
under the two components is caused by 'passage' of the same 
emission entity through our sightline as it rotates around the 
magnetic axis, the average profile-component separation in 
longitude of about 20° suggests a P2 value of about 11° 
for m = 1 — which is consistent with the estimates at meter 
wavelengths. Further, the azimuth angle subtended by the 
profile components at the magnetic axis can be computed 
using our knowledge of the viewing geometry [referring to 
eq. (3) of Paper-I], and is estimated to be about 33°. We 
adopt here the values of a and /3, determined in Paper-I — 
that is, 11.64° and -4.31°, respectively. It follows that the 
magnetic azimuthal separation between the subbeams would 
be about 18° (corresponding to P2), suggesting a rotating 
system of about 20 subbeams. 

It is clear that the 35-MHz observations exhibit phase 
modulation very similar to that noted at higher frequencies, 
and the true modulation frequency is near 0.541 c/Pi, im- 
plying that P3 = 1.848Pi. This means that the subpulse 
drift proceeds in the direction opposite as that of star's ro- 
tation. Some enhancement in the spectral power is evident 
at 0.06 c/Pi, corresponding to the aliased second harmonic 
of the primary modulation feature. The frequency and Q- 
value of the fluctuation features are in excellent agreement 
with those based on analysis of 430-MHz observations. 

The feature 2, at about 0.027 c/Pi, is new and also has 
a high Q. As both LRF- and HRF-spectra show the feature 
at the same frequency, we conclude that it is not aliased 
and that 0.027 c/Pi is its true frequency. It is important to 
note that this feature, seen in set G, is not apparent in the 



corresponding spectra of 111 and 430-MHz data (Paper-I), 
as well as in our other data sets at 35 MHz. We discuss its 
significance below. 

The Amplitude-Modulation Feature, Subpulse 
Spacing and the Number of Rotating Subbeams 

With an understanding of the intrinsic frequencies of the two 
principal features in the LRF spectrum, we can now explore 
a possible relationship between them. A straightforward cal- 
culation, 0.541/0.027 = 20, shows that the two features are 
harmonically related within their errors. We have looked for 
the sidebands (at ± 0.027 c/Pi) around the phase modula- 
tion feature at 0.541 c/Pi. Such side-bands were identified 
in some sections of the 430-MHz sequences (see Paper-I and 
DRa) , but they occur in our data only at a much lower level 
of significance (e.g., see the central panel of Figure 1). The 
0.027 c/Pi feature, just as did the 430-MHz sidebands, sug- 
gests a slower amplitude modulation, with a periodicity, P3, 
of 37.5 Pi or just 20 P3. The magnetic azimuth angle cor- 
responding to the longitude interval Pi would be ~ 18°, if 
a total of 20 subbeams were to be distributed more-or-less 
evenly around the magnetic axis, a value which is in excel- 
lent agreement with our estimate above, implying consis- 
tency with the viewing geometry as determined in Paper-I. 

All these circumstances, along with the high Qs of the 
modulation features and their harmonic relationship, im- 
ply a stable underlying pattern of 20 subbeams, with the 
period of the low frequency feature (feature 2) providing 
a direct estimate of the circulation time of this pattern. 
This slow modulation may also be partially a phase modula- 
tion, considering the asymmetry in the features at 0.027 c/Pi 
and 0.973 c/Pi (see Figure 2). This would imply some non- 
uniformity in the spacing of the subbeams. 

Figures 6 and 7 show the average polar emission pat- 
terns corresponding to the G and C observations, respec- 
tively. These patterns were reconstructed using the 'carto- 
graphic transformation' technique developed in Paper-I as 
well as the various parameters describing our viewing ge- 
ometry of the pulsar determined there, but take all values 
of the circulation time (P3) from the present analysis. In 
assembling these maps, we make no use of the subbeam 
number once the circulation time has been determined di- 
rectly; rather, the 'cartographic transformation' transforms 
the pulse sequence into the frame of the pulsar polar cap 
[see Paper-I : eqs. 5-8], so that the elements of the emit- 
ting pattern and their distribution can be viewed directly. 
In addition, the inverse 'cartographic' transform provides a 
useful means of assessing the consistency and uniqueness of 
the input parameters (see Deshpande, 1999 for details). In 
spite of our overall poor signal-to-noise observations, these 
'closure' consistency checks show a remarkable sensitivity to 
the circulation time, the sign of /3, and to the drift direction. 
Relative insensitivity to the other parameters, however, does 
not allow us to assess the viewing geometry independently. 

Polar Emission Patterns at 35 MHz 

Figure 6 shows the average emission pattern mapped us- 
ing 512 pulses from the G sequence. Remarkably, the 'ac- 
tive' subbeams seem to occupy only about half of the conal 
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Figure 4. Diagrams showing the results of folding the C obser- 
vation at the primary modulation period P3. Note the separate 
'drift' bands for two components in the central panel. The average 
profile is plotted in the bottom panel, and the average flux at a 
particular phase of modulation in the left-hand panel. Note that 
the 'drift' band for the trailing component peaks earlier than that 
for the leading component. 
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Figure 6. Average polar emission map made with 512 pulses 
from the G observation. Note the strong asymmetry about the 
magnetic axis. See text for details. 
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Figure 5. The phase of the primary modulation P3 as a function 
of pulse longitude for the C observation, along with the average 
(solid curve) and modulated-intensity (dashed curve) profiles. 



periphery. Such a deep and systematic modulation in the 
subbeam intensity around the emission cone is completely 
consistent with the appearance of the slow modulation fea- 
ture in the fluctuation spectra (see Figures 1 and 2). The 
emission pattern of the C observation in Figure 7, on the 
other hand, exhibits most of the 20 subbeams and a more- 
or-less uniform spacing — that is, an unsystematic pattern 
of intensities — a picture which is also entirely compatible 
with the observed absence of any significant slow modula- 
tion feature in its fluctuation spectra (Figure 3). It is worth 
mentioning that the power associated with the primary mod- 
ulation (i.e. drift) is significantly higher in the C sequence 
than in the G observations. 

It is clear from both the images that the sampling of the 
subbeams in the radial direction is more or less complete. 
This is significantly so in contrast with the partial sampling 
of the subbeams at 430 MHz. The situation is easily un- 
derstood in terms of the larger angular radius of the radia- 
tion cone at low radio frequencies, which permits our fixed 
sightline to sample it more centrally. The slight elongation 
(in the azimuthal direction) apparent in the subbeam shape 
may not be real, perhaps resulting simply from slight jitter 
in the positions of the subbeams relative to the average. 

A series of polar maps, each from successive subse- 
quences (of duration > P3), were generated and viewed in 
a movie-like fashion^ Although individual frames have a re- 
duced signal-to-noise ratio, certain significant characteristics 
can nonetheless be readily discerned. While some features 
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Figure 7. Average polar emission map made with 140 pulses 
from the C observation. Here we see a nearly full complement of 
subbeams, distributed evenly about the magnetic axis. See text 
for details. 

are stable over only a few circulation times, most {e.g., the 
bright subbeams of the G sequence in Figure 6) remain rel- 
atively prominent over many hundreds of pulse periods. In 
both observations, in a set of successive frames of a movie, a 
few subbeams usually dominate others in intensity. Overall, 
the brightness of the subbeams was observed to fluctuate by 
up to about a factor of 4 over the length of the sequence. 

In the G sequence a few subbeams are persistently 
bright, accounting for most of the power in the average pro- 
file. A very different behavior is seen in the C observation, 
where the subbeams are all very weak and unstable in the 
initial part of the sequence. Then, they gradually increase in 
intensity until most of the 20 subbeams constitute a more- 
or-less stable configuration — that depicted in Figure 7 — for 
some four or five circulation times, whereupon they again 
gradually fade back into a disorganized pattern. The phase 
modulation feature is not readily detectable in either the 
pre- or post-bright phases. Neither, interestingly, are asso- 
ciated changes observed in the shape of the average profile; 
rather, it maintains a clear 'B'-mode profile form, apart from 
the changes in the average intensity. 



3 DISCUSSION 

In the foregoing paragraphs, we have outlined our analysis of 
the fluctuation properties apparent in the decametric emis- 
sion of pulsar B0943+10, as well as our efforts to understand 
the underlying emission pattern responsible for them. 

Although the signal-to-noise ratio of our 35-MHz ob- 
servations is insufficient to detect most individual subpulses, 



we have demonstrated that the average fluctuation behavior 
can nonetheless be reliably characterized, because both the 
magnitude of the fluctuations is large and their pattern rel- 
atively stable. Indeed, such long sequences of pulses present 
an advantage (and opportunity) only if the underlying mod- 
ulation is relatively stable over the relevant time scales, as in 
the present case. The primary phase-modulation frequency 
is observed to vary, but only on scales of many hundreds of 
pulse periods. Although the observed variations in the ap- 
parent P3 might suggest corresponding fractional changes 
in the circulation time, they are more likely to be a mani- 
festation of tiny (an order of magnitude smaller) fractional 
changes in the subbeam spacing. 

Pulsar B0943+10 is unusual in that conal narrowing at 
frequencies above about 300 MHz causes its beam to pro- 
gressively miss our sightline. For a given f3, we expect our 
sightline to sample the radiation cone at different f3/ p values, 
such that more central traverses occur at lower frequencies — 
the 'single' profile form of the pulsar at 430 MHz thus evolv- 
ing smoothly into a well defined 'double' form as f3/p de- 
creases sufficiently at longer wavelengths. The average pro- 
file of PSR 0943+10 at 35 MHz exemplifies this situation 
with an estimated /3/p of 0.78 (see Table II of Paper-I), 
corresponding to a nearly complete sampling of the radial 
'thickness' (Ap) of the hollow-conical emission pattern. The 
fractional thickness of this conal beam Ap/p is estimated to 
be about 15%, consistent with that estimated by Mitra & 
Deshpande (1999). 

The average profile at decameter wavelength also shows 
asymmetry (the peak intensities of the two components 
around the apparent central longitude differ by a factor 
of about 2), so also the profile of modulated power (see 
Figure 5). Such an asymmetry, though striking, may be 
traceable to relatively small departures of the polar cap 
boundary from circularity, if the departures are compara- 
ble to the fractional thickness of the emission cone (Ap/p). 
A study, by Arendt & Eilek (2000), of the detailed shape 
of the polar-cap boundary and its dependence on the in- 
clination angle a, shows that an about 5% departure from 
circularity is expected for a — 15°. Our results, together 
with those of Paper-I, should allow us to estimate the frac- 
tional non-circularity which would cause the observed asym- 
metries, and our efforts in this regard will be reported else- 
where. 

The main results of our decametric study of 
PSR 0943+10's fluctuation properties are as follows: 

• We have shown that fluctuation spectra at decameter 
wavelengths also display a narrow feature with a high Q, 
which is related to the 'drifting' character of its subpulses. 

• We have independently resolved how this feature is 
aliased and have determined the true frequency of the pri- 
mary phase modulation. 

• The high and low frequency features in the G sequence 
were found to be related harmonically (the former just 
20 times the latter). Furthermore, weak sidebands appear 
to be associated with the primary phase-modulation fea- 
ture, again suggesting that a slower fluctuation, of period 
Pz = 2OP3, is modulating the primary phase fluctuation. 

• Since the viewing geometry at 35 MHz is too central 
(J3j p = 0.78) to follow the 'drift' band continuously across 
the pulse window, we have adopted different means of esti- 
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mating the value of P2, which at 11° agrees well with that 
determined in Paper-I. 

• Our observations lead to the conclusion that a system of 
20 subbeams, rotating around the magnetic axis of the star 
with a circulation time P3 of about 37 Pi, is responsible for 
the observed stable subpulse-modulation behavior. 

• Compared with the higher frequency maps, the 35-MHz 
maps sample the subbeams much more completely. However, 
the 35-MHz subbeams, on the whole, show much less uni- 
formity in their positions and intensities. 

• At both 35 MHz and the higher frequencies, the inten- 
sities of individual subbeams fluctuate, maintaining a stable 
brightness only for a few circulation times. The stability 
time-scales of the entire pattern appear roughly comparable 
between 35 MHz and 430 MHz. 

Since the particle plasma responsible for the pulsar 
emission originates close to 'acceleration' zone (which is con- 
sidered to be located near the magnetic polar cap), whereas 
the conversion to radio emission occurs at a height of sev- 
eral hundred kilometers above the stellar surface, we are led 
to conclude that the polar cap is connected to the region 
of emission via a set of 'emission columns'. This picture is 
further supported by the very similar emission patterns ob- 
served at 35 MHz and 430 MHz — emitted, presumably, at 
significantly different altitudes in the star's magnetosphere. 
If indeed some 'seed' activity results in the common charac- 
ter of the emission patterns at different frequencies, then the 
parameters describing the fluctuations (e.g. P2,Pa) should 
be independent of the observing frequency, e.g. as in the 
model of Ruderman & Sutherland (1975). This expectation 
is consistent with our observations of B0943+10, both here 
and in Paper-I, as well as with the findings of Nowakowski 
et al. (1982). 

We expect this basic picture of subbeams and their ap- 
parent circulation to be valid in other pulsars too, although 
the viewing geometry and other quantitative details may dif- 
fer. Whether we observe an amplitude or phase modulation 
would depend largely on the viewing geometry. For exam- 
ple, at frequencies much lower than 35 MHz, we expect the 
phase modulation to be less discernible in B0943+10, and 
the fluctuation would be describable as primarily an ampli- 
tude modulation within each of the subbeam passages across 
our sightline. We suspect that the situation at 35 MHz may 
not be far from such a transition, and therefore any triv- 
ial interpretation of the 'apparent' drift-band separation as 
P2 would be completely misleading. Clearly, a different ap- 
proach is then needed to determine P2 in such a case, as 
was done with our data, or for pulsar B0834+06 (Asgekar 
& Deshpande, 2000). 

Among a total of 6 detections at 35 MHz, we find only 
one sequence exhibiting a 'Q'-mode behavior — that is, with 
both the average profile and fluctuation spectra displaying 
clear deviations from 'B'-mode properties. This predomi- 
nant 'B'-mode character of the decametric emission in our 
observations, as well as in the earlier observations of Desh- 
pande & Radhakrishnan (1994), appears inconsistent with 
the trend at higher frequencies where both the 'B' and 'Q' 
modes may occur with similar frequency (Suleymanova et 
al. 1998). It is possible however, that the intervals of dis- 
organized pattern in our observations correspond to some 
mode-change activity. The weakening of the pulse after the 



bright phase in the C sequence may be regarded as such an 
example. 

To summarize, our observations of PSR 0943+10 at 
35 MHz show remarkable similarity with the high fre- 
quency observations in Paper-I in terms both of their fluc- 
tuation properties and of the underlying emission patterns 
they represent. This similarity, combined with the 'radius- 
to- frequency mapping', implies that the regions of radio 
emission correspond to a well-organized system of plasma 
columns in apparent circulation around the magnetic axis of 
star. We emphasis that the observed circulation should not 
be mistaken for an unlikely physical (transverse) motion of 
the plasma columns in the presence of strong magnetic fields. 
The circumstances suggest that a common 'seed' activity is 
responsible for the generation and motion of the relativistic 
plasma. 
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